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ABSTRACT 

Context. The dusty nuclear regions of luminous infra-red galaxies (LIRGs) are heated by either an intense burst of massive star 
formation, an active galactic nucleus (AGN), or a combination of both. Disentangling the contribution of each of those putative dust- 
heating agents is a challenging task, and direct imaging of the innermost few pc can only be accomplished at radio wavelengths, using 
very high- angular resolution observations. 

Aims. We observed the nucleus A of the interacting starburst galaxy Arp 299, using very long baseline interferometry (VLBI) radio 
observations at 1.7 and 5.0 GHz. Our aim was to characterize the compact sources in the innermost few pc region of Arp 299-A, as 
well as to detect recently exploded core-collapse supemovae. 

Methods. We used the European VLBI Network (EVN) to image the 1.7 and 5.0 GHz compact radio emission of the parsec-scale 
structure in the nucleus of Arp 299-A with milliarcsecond resolution. 

Results. Our EVN observations show that one of the compact VLBI sources, Al, previously detected at 5.0 GHz, has a flat spectrum 
between 1.7 and 5.0 GHz and is the brightest source at both frequencies. Our 1.7 GHz EVN image shows also diffuse, low-surface 
brightness emission extending westwards from Al and displays a prominent core-jet structure. 

Conclusions. The morphology, radio luminosity, spectral index and ratio of radio- to-X-ray emission of the A1-A5 region is consistent 
with a low-luminosity AGN (LLAGN), and rules out the possibility that it is a chain of young radio supemovae (RSNe) and supernova 
remnants (SNRs). We therefore conclude that A1-A5 is the long-sought AGN in Arp 299-A. This finding may suggest that both 
starburst and AGN are frequently associated phenomena in mergers. 
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1. Introduction 

Arp 299 consists of two interacting galaxie s (IC 694 and NGC 
3690) that are in an early merger sta ge (|Keel & Wu| [1995). 
At a luminosity distance of 44.8 Mpc dFixsen et al.l 119960 for 
\ H = 73 km s _I Mpc" 1 (1" corresponds to 217 pc) , Arp 299 ha s 
an infrared luminosity, Lir as 6.7 x 10 n L o dSanders et alJl2003l) . 
The nuclear regions of Arp 299 are heavily dust-enshrouded, 
which prevents their study at optical wavelengths. Fortunately, 
radio observations are not affected by dust, and therefore stand 
as the most promising way to study the innermost, buried regions 
of Arp 299. 

Radio interferometric VL A observations o f Arp 299 at about 
one arcsecond resolution (e.g. Neff et al. 2004) show that its ra- 
dio emission is mainly due to ionized gas. The large distance to 
Arp 299 requires the use of very long-baseline interferometry 
(VLBI) technique to resolve the radio emission from the inner- 
most regions of this interacting starburst galaxy. In fact, VLBI 
observations of the two primary nuclei of Arp 299 have resulted 
in the discovery of a large population of very compact sources 
in the nuclear regi ons of Arp 299-A (-IC694) and Arp 299 - 
Bl (=NGC 3690) dPerez-Torres et alJl2009bt Iuivestadll2009h . 
The high brightness temperatures observed for this population of 
compact sources were indicative of a nonthermal origin for the 
observed radio emission, implying that most of those sources 
were young RSNe and SNRs . On the other hand, X-ray imag- 
ing and spectrosc opy dDella Ceca et al. 20021: IZezas et alj f2003; 
iBallo et all2004T) . as well as the presence of nuclear H2O masers 



(iHenkel et alJl2005l:lTarchi et al.ll2007l) . indicate that the nuclear 
regions of Arp 299-A and Arp 299-B 1 are likely to contain an 
AGN, but its lo cation has so far remained an open question. On 
the other hand, lAlonso-Herrero et alJ (l2009h have noticed that 
the MIR properties of Arp 299-A were consistent with obscured 
star formation, although they could not exclude the existence of 
an LLAGN. In this Letter, we focus on studying the inner few 
parsecs of the nucleus in Arp 299-A. A detailed discussion of the 
radio emission and spectral properties of all the compact com- 
ponents in the central hundred parsec region will be presented in 
another paper. 

2. EVN observations of Arp 299-A 

We used the European VLBI Network (EVN) to image the cen- 
tral region of Arp 299-A at 1.7 GHz on 7-8 June 2009, and at 
5.0 GHz on 12-13 June 2009, with the pointing position centered 
on the Arp 299-A nucleus. We used a sustained data recording 
rate of 1024 Mbit s in two-bit sampling. Each frequency band 
was split into eight intermediate frequencies (IFs) of 16 MHz 
bandwidths each, for a total synthesized bandwidth of 128 MHz. 
Each IF was in turn split into 32 channels of 0.5 MHz bandwidth 
each. Our array included the following antennas at both fre- 
quencies: Effelsberg, Cambridge, Medicina, Jodrell Bank (Mk 
2), Onsala, Torun, Urumqi, Shanghai, and Westerbork. In addi- 
tion, the Yebes antenna in Spain was also part of the array at 5.0 
GHz. Each of the two VLBI observing epochs consisted of 6.0 hr 
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phase-referenced experiments. Our target source, the nucleus of 
Arp 299-A, was phase-referenced to the calibrator Jl 128+5925 
with a duty cycle of 5 minutes. The strong source 4C 39.25 was 
used as fringe finder and bandpass calibrator. The total effective 
time on source for Arp 299-A was ~3.6 hr. 

The data were correlated at the EVN MklV data processor 
at JIVE using an averaging time of 2 s. We performed stan- 
dard a-priori gain calibration using the measured gains and sys- 
tem temperatures of each antenna. This calibration, as well as 
the data inspection and flagging, were done within the NRAO 
Astronomical Image Processing System (AIPS). We also cor- 
rected for ionosphere effects (of particular relevance at 1 .7 GHz) 
and source-structure effects of the phase-reference source at 
both frequencies, followi ng the same procedures as described in 
Pere z-Torres et al. (2009b|). No self-calibration was performed 
on the data, since the peaks of emission were too faint to qualify 
for these procedures. We also used AIPS to image Arp 299-A, 
that is, the images in Figure [TJ were obtained by applying natu- 
ral weighting to the data, which resulted in FWHM, synthesized 
interferometric beams of (11.9 mas x 4.4 mas at PA=-76°) and 
(5.0 mas x 4.1 mas at PA=7°) at 1.7 and 5.0 GHz, respectively. 
The attained \<j off-source r.m.s noise was of ~ 25yt/Jy at 1.7 
GHz and ~ 22 /Jy at 5 .0 GHz. 

3. Results and discussion 

Our results are summarized in Figure [TJ and Table QJ Figure [1] 
shows the field of the central (150x80) pc of Arp 299-A (top), 
as imaged with the EVN at 5.0 GHz on 12-13 June 2009. Our 
5.0 GHz E VN image is very s i milar t o the images previously 
published in IPerez-Torres et al.l (l2009bl) . and confirms that Arp 
299-A hosts a larg e number of bright , com pact radio sources. 
As already noted bv lPerez-Torres et all d2009bl) . the high bright- 
ness temperatures displayed by those compact components indi- 
cate that they must be nonthermal emitting sources, most likely 
young radio SNe and SNRs, and a detailed discussion of the ra- 
dio emission and spectral properties of all these compact sources 
will be presented elsewhere. The middle and bottom panels of 
the figure are blow-ups of the inner 8 parsec nuclear region of 
Arp 299-A, as imaged with the full EVN at 1.7 and 5.0 GHz, 
and the positions of its compact components, as well as their 
flux densities, are summarized in Table QJ 

The main results from our observations concern the nature of 
component AO and, especially, of the complex formed by A1-A5 
(lower panel in Figure [1] and Table QJi, which we discuss in the 
next sections. 

3. 1 . The young radio supernova AO 

Component AO was identified with an RSN bv lNeff et all (l2004h . 
AO showed a strongly inverted spectral index between 2.3 and 
8.4 GHz, indicating that it was a very young SN. We al- 
ready detected it at 5.0 GHz dPerez-Torres et alll2009bl) . more 
than five years after its discovery at 8.4 GHz, which implied 
that AO was a longlasting , slowly evolving, nonstandard RSN 
(IPerez-Torres et al.ll2009M) . Our current 5.0 GHz EVN obser- 
vations show that AO continues to be very bright, with L5.0 = 
7.4x 10 26 erg s _1 Hz -1 , confirming that it is a longlasting, slowly 
evolving Type II supernova in a rich, dense surrounding medium. 
Not detecting AO at 1.7 G Hz also confirms the suggestion by 
IPerez-Torres et alJ (l2009bl) that its radio emission is being sup- 
pressed at low frequencies by a nearby absorber, most likely 
a foreground H II region in the vicinity of AO. We can char- 
acterize the main parameter of this H II region, its emission 



measure EM, by assuming an (intrinsic) spectral index typi- 
cal of Type II SNe in their optically thin phase [a = -0.75; 
5mtr,i.7 = 55.o(l-7/5.0) u ']; and purely foreground, free-free ab- 
sorption (S bs,\.i = 5intr,i.7 e~ TBX1 ). We then obtain Tffj.7 > 2.2. 
The EM at GHz frequencies is proportional to the free-free opti- 



cal depth and is given by EM = 3.1x10 r ff v T 



1.35 ,,2.1 



cm pc 



dMezger & Hendersonll9 67). where T e ^ is in units of 10 4 K, and 
v in GHz. Assuming the free-free absorber is a purely hydrogen 
plasma with T e ^ = 1, we obtain EM > 2.0 x 10 7 crrT 6 pc for 
the H II close to AO. This value of EM is very similar to the one 
obtained fo r the foreground H II region close to the well studied 
SN 2000ft dAlberdi et alj2006l;lPerez-Torres et al.ll2009al) . in the 
dust-enshrouded galaxy NGC 7469. 



3.2. The A1 -A5 complex 

We detected compo nents Al and A5, which w ere already re- 
ported at 5.0 GHz bv lPerez-Torres et al.ld2009bl) . Al, the bright- 
est compact source detected during our previous 5.0 GHz eEVN 
observations, continues to be the brightest one at 5.0 GHz 
(L5.0 = 1.8 x 10 27 erg s _I Hz -1 ), and is also detected at 1.7 GHz 
{L\ j = 2.0 x 10 27 er g s' 1 Hz" 1 ). Component A5 is also detected 
at both frequencies. IPerez-Torres et al . (2009b) suggested that 
Al could be an SNR, while A5 it was suggested a possible SN, 
based on its marginal variability at 5.0 GHz (although we could 
not exclude that A5 was an SNR). However, our new, contempo- 
raneous observations at 1.7 and 5.0 GHz shed new light on their 
nature, and suggest a rather different scenario. 

In fact, Al and A5 are part of a "global structure", as clearly 
shown by our 1.7 GHz EVN observations (see lower panel in 
FigureQJi, where a bridge of radio emission connects them, with 
two new regions of enhanced emission in between (components 
Alb and Ale) . The extension of the A1-A5 complex is ^ 31x21 
mas, corresponding to a linear size of ^ 6.7x4.6 pc. Within the 
1.7 GHz structure, we detect two new components, Alb and 
Ale, which are not detected at 5.0 GHz more than three times 
the off-source r.m.s. Nevertheless, detecting them clearly at 1.7 
GHz (> 9cr) implies that they are real, steep-spectrum compo- 
nents, whose nature we discuss in the next two sections. 



3.2.1 . A chain of nested SNe and SNRs in an SSC...? 

The spectrum of Al between 1.7 and 5.0 GHz is very flat 
(a = -0.13 + 0.11). This spectral index would in principle be 
compatible with an H II region but, as we already showed in 
Pere z-Torres et al.|d2009bl) . all compact components in Arp 299- 
A show brightness temperatures that imply a nonthermal origin 
for their radio emission. Thus, the flat spectrum shown by Al 
could be explained by radio emission from an SNR or, alterna- 
tively, from an AGN. How can we distinguish these two scenar- 
ios? To answer this question we must consider Al in the context 
of the complex it is part of, together with components Alb, Ale, 
and A5. As with Al, the brightness temperatures of the other 
components rule out a thermal origin for their radio emission. 
In principle, the spectral indices of those components (see Table 
[TJ are compatible with Alb and Ale being young supernovae 
in their optically thin phase, and A5 being another SNR, all of 
them part of a super star cluster (SSC). 

The existence of SSCs in Arp 299-A has been demon- 
strate d by their det ection using 2.2-/im adaptive optics imag- 
ing dLai et alJ [l999). and further evidence comes from Hubble 
Space Telescope (HST) FOC and NICMOS images, which reveal 
a population of young (4-15 Myr) stellar clusters in the cen- 
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Source 


Source 


Aa" 


AS* 


Angular size b 


S1.7 G"Jy) c 


s 50 (^Jy) c 


Spectral index 


Name 


Type 


(J2000.0) 


(J2000.0) 


(milliarcsec) 






a 50 

"l 7 


AO 


SN 


33.6213 


46.707 


< 4.1 x 4.0 


<75 


307±27 


>1.28 


Al 


AGN (core) 


33.6199 


46.699 


< 6.1 x 2.7 


822±48 


716±42 


-0.13+0.11 


Alb 


AGN (jet) 


33.6190 


46.699 


< 10.1 x 5.1 


314±30 


<66 


<-1.42 


Ale 


AGN (jet) 


33.6183 


46.695 


< 14.3 x 6.2 


267±28 


<66 


<-1.27 


A5 


AGN (jet) 


33.6176 


46.694 


< 6.9 x 5.0 


246±28 


183±24 


-0.27±0.22 



a Coordinates are given with respect to a= 1 1:28:00.0000 and 6 = 58:33:00.000, and the typical uncertainty is less than 1 mas. 

b Upper limits to the deconvolved angular size (Major axis x Minor axis, in mas) as obtained using task IMFIT within AIPS, using our 5.0 GHz 
for all components but Alb and Ale, which were only detected at 1.7 GHz. 

c The uncertainty in the reported flux density (and derived spectral index) for the detected compact components corresponds to lcr, where cr was 
determined by adding in quadrature the off-source r.m.s. in each image and a 5% of the local maxima, to conservatively account for possible 
inaccuracies in the EVN calibration. Upper limits for undetected sources are three times the off-source r.m.s. 



tral regions of Arp 299-A dAlonso-Herrero et al J [20001) . High- 
resolution near-infrare d spectra of the nuclear starburst of M 82 
dMcCradv et al.ll2003l) have shown that its SSCs have typical 
radii of 2-3 pc, and masses that range from ~ 10 5 M Q up to (1- 
2)xl0 6 M Q . The size of the A1-A5 complex is 6.7 x 4.6 pc, in 
agreement with expectations for SSC sizes. Therefore, a possible 
scenario is one where a young SSC has resulted in two young, 
recently exploded SNe (Alb and Ale), and two SNRs (Al and 
A5). In this case, the steepness of the spectral indices for sources 
Alb and Ale can only be explained as from very recently ex- 
ploded (f < 10-20 yr) core-collapse supernovae. Furthermore, 
to estimate the age of the SNR candidates Al and A5, we use up- 
per limits to their (deconvolved) angular size (column 5 in Table 
[TJ. Al and A5 have a maximum angular radius of 6.1 and 6.9 
mas, which correspond to a linear size of 1.3 and 1.5 pc, respec- 
tively. Assuming an average expansion speed of 3000 km s~'for 
both SNRs, those sizes correspond to kinematic ages of ^430 yr 
for Al and ^490 yr for A5. 

In summary, if the A1-A5 complex is an SSC, it has yielded 
four bright radio supernovae in the past 500 yr. Now, this 
is very unlikely to happen in an SSC of just 3.4 pc radius 
whose mass is unlikely to be higher than (l-2)xl0 6 M . Indeed, 
evolutionary models for th e radio emission in starbursts by 
iPerez-Olea & Colinal d 19951) show that the SN rate in an instan- 
taneous burst can be as high as 20xl0~ 10 yr _1 M _1 up to an age 
of ^6 Myr, and about 10xl0~ 10 yr~ lM o~' at an a g e of 9-15 Myr. 
Thus, an SSC of 1O 6 M would yield <1 SN every 500 yr at an 
age of 6 Myr, and <1 SN every 1000 yr at an SSC age of 9-15 
Myr. Therefore, if A1-A5 were SNe and SNRs exploding in a 
SSC, we should have expected to detect at most one SN/SNR in 
the past 500 yr, whereas we find four. The probability of detect- 
ing two (young) SNe in the past 20 yr, and two SNRs with ages 
of -430 and ~ 490 yr in the past 500 yr is less than 3 x 10~ 6 . All 
this evidence seems to rule out an SSC scenario for the A1-A5 
complex. 



3.2.2. ...Or the long-sought AGN in Arp 299-A? 

We show in the next paragraph that a much more plausible sce- 
nario for the A1-A5 complex is that of an AGN, where Al-the 
brightest compact source at both frequencies in the nuclear and 
circumnuclear regions of Arp 299-A- is identified with the core 
of the AGN in Arp 299-A, while components Alb, Ale, and A5 
are jet components (see Figure [TJ. There are strong grounds for 
this scenario: 



(i) A1-A5 displays a core-jet morphology, typical of AGNs. 

(ii) The two-point spectral indices of the A1-A5 complex (see 
Table [TJ agree well with expectations from an AGN. In this 
scenario Al, which has the flattest spectrum between 1.7 and 
5.0 GHz, would be the core, and Alb, Ale, and A5 would be 
jet components, (iii) There is evidence of strong flux density 
variability of A 1 at 2.3 and 8.4 GHz betwee n 2002 and 2003 
dNeff etalJl200l IUlvestadll2009l) . Figure 5 of lUlvestadl (2009) 
shows little variation for Al starting from 2004, or even a slight 
decrease (increase) in flux density at 2.3 (8.4) GHz. On the 
other hand, there is no evidence for significa nt variability at 5.0 
GHz b etween April 2008 and June 2009 (see Perez-Torres et al.l 
2009b and Table [TJ. This behavior is at odds with an SN ori- 
gin for the radio emission of Al, and is more easily reconciled 
withi n an AGN scen ario. In addition, by combining the data 
from Ulvestad (2009) with ours and assuming no variability be- 
tween 2005 and 2009 at cm-wavelengths, we find that the spec- 
trum of Al must have a turnover frequency between 1.7 and 5.0 
GHz, as expected for the core of an AGN whose radio emission 
is partially self absorbed, (iv) The radio luminosity of the A1-A5 
complex between 1.7 and 8.4 GHz is ^ 1.9 x 10 37 erg s _1 , typ- 
ical of LLA GNs. (v) The ratio of (5 G Hz) radio-to-X ray (soft) 
luminosity (Terashima & Wilson 2003) for the A1-A5 complex 
in Arp 299-A is typical of LLAGNs. From Table [TJ we ob- 
tain vL Vj 5.o ~ 2.0 x 10 37 ergs -1 . The combined absorption 
co rrected X-ra y lumi nosity for the X-ray sources discovered 
bv IZezas et aT] d2003l) in Arp 299-A is ~ 1.9 x 10 4() erg s _1 . 
This results in a value of Rx — 1.1 x 10~ 3 , which agrees 
closely with expectations from an LLAGN (see Figure 4 of 
iTerashima &~W ilson 2003|). 



We therefore identify the A1-A5 complex with the AGN in 
Arp 299-A. In this case, the cm-radio luminosity of the AGN in 
Arp 299-A is very similar to the one displayed by the LLAGN 
in M 81, which is the closest galaxy with an active nucleus. 
The core of M 81 also presents a faint, one-sided jet, with a 
flat spectrum near the core and a steep spectrum along the jet 
(Biet enholz et al.l2004l) . Furth ermore, cm-VLBI observatio ns of 
the merger system NGC 6240 dGallimore & Beswickll2004l) and 
mm-Plateau-de-Bure inter ferometer observations o f the proto- 
typical ULIRG Arp 220 dDownes & Eckartl l2007h have sug- 
gested the coexistence of both AGN and starburst in their nu- 
clear regions. Thus, our discovery suggests that both starburst 
and AGN are frequently associated in merging systems. 
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4. Summary 

We imaged with milliarcsecond resolution the inner eight par- 
sees of the nuclear region in Arp 299-A, using contemporane- 
ous EVN observations at 1.7 and 5.0 GHz. At 5.0 GHz, we 
detect components AO, A l, and A5, as previously detected by 
iPerez-Torres et al.1 (l2009bl) . At 1 .7 GHz, Al and A5 are also de- 
tected, but not AO. In addition, two new components between 
Al and A5 (Alb and Ale) are detected at 1.7 GHz, forming a 
complex. 

We find that the morphology, spectral index, radio luminos- 
ity, and radio-to-X ray luminosity ratio of the A1-A5 complex 
are consistent with that of an LLAGN, and rules out the possi- 
bility that it consists of a chain of young RSNe and SNRs in a 
young SSC. We therefore conclude that Al is the long-sought 
AGN in Arp 299-A. Since Arp 299-A had long been thought of 
as a pure starburst, our finding of a buried, low-luminosity AGN 
in its central region, coexisting with a recent burst of starforma- 
tion, suggests that both a starburst and AGN are frequently as- 
sociated phenomena in mergers. In this case, our result is likely 
to have an impact on evolutionary scenarios proposed for AGN 
and the triggering mechanism of activity in general. 

Finally, we also note that component AO, previously identi- 
fied as a young RSN, is not seen at our low-frequency observa- 
tions, which implies there is a foreground absorbing H II region. 
It is remarkable that this RSN exploded at the mere distance 
(projected) of two parsecs from the putative AGN in Arp 299-A, 
which makes this supernova one of the closest to a central super- 
massive black hole ever detected. This result may also be rele- 
vant to accreting models in the central regions of galaxies, since 
it is not easy to explain the existence of very massive, supernova 
progenitor stars so close to an AGN. While seemingly contradic- 
tory, this could explain the low-luminosity of the AGN we see 
in Arp 299-A. In fact, since massive stars shed large amounts of 
mechanical energy into their surrounding medium, thereby sig- 
nificantly increasing its temperature, those massive stars would 
hinder the accretion of material to the central black hole, which 
could in turn result in a less powerful AGN than usual. 
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Fig. 1. Top: 5.0 GHz full EVN image of the central 150 parsec 
region of the luminous infrared galaxy Arp 299-A (=IC 694), 
displaying a large number of bright, compact, nonthermal emit- 
ting sources, mostly identified with young RSNe and SNRs. The 
image center is atRA 1 1:28:33.63686 and DEC 58:33:46.5806. 
Middle and bottom: Blow-ups of the inner 8 parsec of the nu- 
clear region of Arp 299-A, as imaged with the full EVN at 1.7 
and 5.0 GHz. The image center is at RA 11:28:33.61984 and 
DEC 58:33:46.7006 in both panels. The morphology, spectral 
index and luminosity of the A1-A5 region are very suggestive of 
a core-jet structure. The color scale goes from -50 /iJy/b up to 
400 /iJy/b in the top panel and from 125 /iJy/b to 400 /uJy/b in 
the middle and bottom panels. Contours are drawn at 5 and 10 
times the off-source r.m.s. noise. 



